Mon. Not. R. Astron. Soc. 000, [[ju] (1994) Printed 1 February 2008 (MN M^jX style file vl.4) 



O 
O 

(N 



Real-Time Difference Imaging Analysis of MOA Galactic 
Bulge Observations During 2000 

LA. Bond, 1 ' 2 F. Abe, 3 R.J. Dodd, 1 ' 4 ' 5 J.B. Hearnshaw, 2 M. Honda, 6 
J. Jugaku, 7 P.M. Kilmartin, 1 ' 2 A. Maries, 1 K. Masuda, 3 Y. Matsubara, 3 
Y. Muraki, 3 T. Nakamura, 8 G. Nankivell, 5 S. Noda, 3 C. Noguchi, 3 
K. Ohnishi, 9 N.J. Rattenbury, 1 M. Reid, 4 To. Saito, 10 H. Sato, 8 
M. Sekiguchi, 6 J. Skuljan, 2 D.J. Sullivan, 4 T. Sumi, 3 M. Takeuti, 11 
Y. Watase, 12 S. Wilkinson, 4 R. Yamada, 3 T. Yanagisawa, 3 
and P.C.M. Yock 1 

^ Wn.r11lt.11 nf Rr'i en re TJn.ilirrtii.tii nt Aiirkln/nd A lurhlnm A New 'Zm.l.n.'n.d. 



^Faculty of Science, University of Auckland, Auckland, New Zealand 
2 Department of Physics and Astronomy, University of Canterbury, Christchurch, New Zealand 
(N ' 3 Solar- Terrestrial Environment Laboratory, Nagoya University, Nagoya 464, Japan 

^School of Chemical and Physical Sciences, Victoria University, Wellington, New Zealand 
■ 5 Carter Observatory, P. O. Box 2909, Wellington, New Zealand 
° Institute of Cosmic Ray Research, University of Tokyo, Tanashi, Tokyo 188, Japan 
^ ''Research Institute of Civilization, Tama 206, Japan 

, 8 Department of Physics, Kyoto University, Kyoto 606, Japan 
9 Nagano National College of Technology, Japan 
j- —i i 10 Tokyo Metropolitan College of Aeronautics, Tokyo 116, Japan 
Tohoku University, Sendai, Japan 
, 12 KEK Laboratory, Tsukuba 305, Japan 



^ C"| Accepted . Received ; in original form 

Or 



o 

c3 



ABSTRACT 

We describe observations carried out by the MOA group of the Galactic Bulge during 
2000 that were designed to detect efficiently gravitational microlensing of faint stars in 
which the magnification is high and/or of short duration. These events are particularly 
useful for studies of extra-solar planets and faint stars. Approximately 17 deg 2 were 
monitored at a sampling rate of up to 6 times per night. The images were analysed in 
real-time using a difference imaging technique. Twenty microlensing candidates were 
detected, of which 8 were alerted to the microlensing community whilst in progress. 
Approximately half of the candidates had high magnifications (Jj 10), at least one 
had very high magnification (<; 50), and one exhibited a clear parallax effect. The 
details of these events are reported here, together with details of the on-line difference 
imaging technique. Some nova-like events were also observed and these are described, 
together with one asteroid. 

Key words: Gravitational lensing: microlensing — Techniques: image processing 



1 INTRODUCTION 

The observation of gravitational microlensing events is be- 
coming established as an important tool in several astro- 
physical endeavours. These include the search for extra- 
solar planets, studies of dark matter, studies of stellar at- 
mospheres, and studies of galactic structure. More than 500 
microlensing events, most of them alerted in real time, have 
been identified by the micro lensing survey grou ps MACHO 
(|Alcocket al. 1993|), OGLE flUdalski et al. 1992|), and EROS 



(Aubourg et al. 199S ) . These groups monitored some mil- 
lions of stars within their fields of view to find the microlens- 
ing events. The photometry measure ments were based on 
profile fitting software such as Dophot ( Schechter et al. 199£ ) 
which determines fluxes by integrating over the profile of a 
given star on a given image. 



While profile fitting techniques have proved to be highly 
successful, the limitations of this type of analysis have been 
recognized for some time. The most serious limitation is the 
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contamination of flux measurements by the "blending" of 
the profiles of neighbouring stars. This problem becomes 
more serious when one encounters crowded regions of stars 
within the fields of view. Blending causes systematic errors 
in phot ometry measurements and subsequent systematic er- 



in New Zealand. The telescope was modified by the MOA 
group to provide a 1?3 field of view at f/6.25 and to en- 
able computer controlled tracking (Abe et al. 1997). The 
data were taken using the wide field camera, MOA-cam2, 



that was constructed by the MOA group (Yanagisawa et al 



rors in the determination of microlensing parameters such 200C). MOA-cam2 is a mosaic camera comprised of three 



as event duration and peak amplification ( Han 1998] ). Alard~ 
(2001) has shown that, using an image subtraction or differ- 



ence imaging analysis (hereafter DIA) to allow for the effects 
of blending, an unbiased reconstruction of the mass function 
in the Galactic bulge should be possible using microlensing 
data alone. 

Another problem with profile fitting analyses is that 
only those microlensing events whose source stars are visible 
and resolvable at baseline are detected. Microlensing events 
that rise from below the observational threshold to be vis- 
ible during times of high magnification are missed. This is 
unfortunate, because these events with faint source stars are 
potentially capable of providing valuable information. The 
mere detection of faint stars in the Galactic Bulge that are 
normally undetectable can clearly provide useful informa- 
tion. For example, if the magnification is high, spectroscopic 
information may also be obtained that is not available oth- 
erwise (Lennon et al. 1996, 1997; Minniti et al.1998). Kane 
and Sahu (2001) recently made a kinematic study of the far 
side of the Galactic bulge using spectroscopy of microlensed 
sources. 

The peaks of high magnification events provide valu- 
able information on extra-solar planetary systems. Griest 
and Safizadeh (1998) showed that the probability for detect- 
ing Jupiter-like planets is nearly 100% for events with peak 
amplifications J> 100 provided good photometric sampling is 
carried out around the times of the peaks, and that the prob- 
ability for detecting Earth-like planets is substantial. They 
also showed the probability for detecting Jupiter-like plan- 
ets is substantial in events with peak amplifications <; 10. 
These predictions have been borne out by observation (Rhie 
et al. 2000; Albrow et al. 2000; Bond et al. 2001). As will 
be shown here, high magnification events can be detected 
efficiently using DIA. 

DIA largely overcomes the problems caused by blended 
or undetected microlensing source stars. Recent re-analyses 
of the MACHO database using DIA have increased the mi- 
crolensing detection rate and improved the quality of the 
photometry (Alcock et al. 1999a; 1999b; 2000). A re-analysis 
of the OG LE database i s underway and similar results are 
emerging ( Wbzniak 200C ) . All of these analyses were carried 
out off-line. 

During 2000 we applied DIA to observations of the 
Galactic Bulge carried out by the MOA Group. The analy- 
sis was "real-time" in that alerts of microlensing events in 
progress were provided to the microlensing community. In 
this paper we describe our implementation of the real-time 
DIA method and present the transient phenomena observed 
during 2000. 



abutted thinned 2kx4k SITe CCD chips. Two wide pass- 
bands were employed: a "blue" (400-630 nm) and a "red" 
(630-1000 nm) band. The wide passbands compensate for 
the small telescope aperture. The total field of view over 
the three CCDs is 55' x 83' and the pixel size on the sky 
is 0.81". This is matched to the less than ideal seeing of 
the set-up which ranges from about 1.8" to 3.5", with the 
median at about 2.5". These are global figures that include 
dome seeing, etc. 

The observational strategy employed by MOA for the 
Galactic Bulge in 2000 was somewhat different from that of 
the established survey groups MACHO, OGLE, and EROS. 
Rather than sample a large area of the Bulge once per night, 
we chose to sample a smaller area several times per night. 
This decision was taken mainly with the aim of improving 
our sensitivity to very high magnification events which are 
expected to rise and fall very rapidly and only be detected 
above the observation threshold for a few days or less. 

For the 2000 observations, a total of 14 MOA- 
cam2 fields towards the Galactic Bulge were selected on the 
basis of the distribution of microlensing events detected by 
the MACHO Collaboration during 1995-1999. Q These fields 
cover an area of ~ 17 deg 2 . The exposure time was 180 sec- 
onds per field with most exposures being taken in the red 
passband. Exposures in the blue passband were only made 
occasionally and are not used in this study. This observa- 
tional strategy allowed a sampling rate for the full 14 fields 
of up to 6 times per night. These fields were observed by 
MOA between 2000 March and September. Problems with 
the camera electronics prevented observations over a six- 
week period during July- August. 



3 DATA ANALYSIS 

3.1 Image Subtraction Method 

Two geometrically aligned images taken under conditions of 
different seeing are related through the convolution relation 



i{x,y) = r * k(x,y) + b(x,y). 



(1) 



Here r is the better seeing image referred to as the "refer- 
ence" image and i is the current "observation" image. The 
convolution kernel k encodes the seeing differences and b 
represents the sky background differences between the two 
images. Solving for the convolution kernel is the crucial step 
in image subtraction. There are two approaches in use for 
doing this. One matches point spread function (PSF) mod- 
els on the two images by solving for the kernel in Fourier 
space (Phillips & Davis 1995; Tomaney & Crotts 1996). This 



2 OBSERVATIONS BY MOA OF THE 
GALACTIC BULGE 

The observations were carried out using the 60-cm Boiler 
& Chivens telescope of the Mt John University Observatory 



* These observation fields are a revision of previously defined 
survey fields observed with MOA-cam2 during 1998-99. The data 
from observation fields used prior to 2000 were not used in the 
survey analysis in the present study, but they were used for follow- 
up analysis of selected events that are presented below. 
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Figure 1. Difference imaging in action. The top left image is a sub-region of a reference image. The top right image is a corresponding 
observation image. The squares in the two top images mark the location of stamps used to calculate the kernel. The lower left image is 
the result of subtracting the reference image from the observation image. The lower right image is the same subtracted image with the 
positions of saturated stars masked out. 



approach was adopted by Alcock et al. (1999a) in their re- 
analysis of a subset of the MACHO data. 

The ot her approach is to dire ctly model the kernel in 
real space (Alard & Lupton 1998). We developed our own 



implementation of this method. This includes the modifica- 



tion of Alard (1999) that models spatial variations of the 
kernel across the CCD. The math emat ical techniq ues ar e 
described fully by Alard & Lupton (tl998ft and Alard (h 999|) . 



Here we summarize our implementation. The convolution 
kernel function at a given location (xo,yo) is expressed as a 
linear combination of analytical "basis" functions 

k(x ,y ,x,y) = a f (x,y) + ^^a m (x ,y )f m (x,y). (2) 



The coefficients in the linear combination are spatially de- 
pendent. They are also modeled as a linear combination of 
basis functions which takes the form of a two dimensional 
polynomial 



i{xo,Vo) 



(3) 



where the polynomial indices p and q depend implicitly on 
the function index n. The basis functions used to model the 
kernel take the form of a combination of two dimensional 
Gaussian functions and polynomials. The function fm(x,y) 
in Eqn. 2 can be expressed in such a way as to render the in- 
tegral of each but the first basis function zero. Consequently 
the integral of the kernel itself is constant across the chip 
regardless of the spatial variations of the kernel. The differ- 
ential background is modeled in the same way as the spatial 
variations in the kernel: 



b(x,y) 



E 



b n x p y 9 



(4) 



where again the polynomial indices depend implicitly on the 
function index. 

The solution for the coefficients a mn and b„ may be 
found using standard linear techniques. The linear system 
is constructed using a number of sub-regions, referred to 
as "stamps", centred on bright, but not saturated, stars 
throughout the image. An advantage of this method of so- 
lution is that these stamps need not be centred on isolated 
stars with clear point spread functions (PSF stars). There 
can be any number of stars to any degree of crowding within 
a stamp. However, we impose the strict requirement that 
regions enclosed by the stamps, plus an additional margin 
given by half the kernel sizeQ, are completely free of bad 
pixels including dead columns and bleeding columns due to 
nearby saturated stars. Imposing this condition greatly im- 
proves the robustness in obtaining good quality subtracted 
images. Another consideration is that the size of the linear 
system is large when one includes the spatial variations in 



t We used a kernel of dimensions 17 X 17 pixels in this analysis. 
© 1994 RAS, MNRAS 000, 



the convolution kernel. The computing time required to cal- 
culate the linear system can be greatly reduced using an 
accelerated summation procedure which ass umes that the 
kernel is constant within individual stamps (Alard 1999). 



Our implementation of the kernel solution method along 
with the subsequent image subtraction is coded in C++. It is 
worth mentioning here the merits of using an object oriented 
programming language for this purpose. Eqns. 2-4 are some- 
what unwieldy to encode in a procedurally oriented language 
such as Fortran or C. The object oriented feature of C++ al- 
lows one to encapsulate neatly a set of basis functions and 
their associated polynomial index tables in a single class. 
One then declares separate class instances for each set of 
basis functions describing the kernel, its spatial variations 
and the differential background. 

The PSF profiles are known to vary across the field for 
large CCD images. Consequently, the possible spatial varia- 
tion of the convolution kernel is an important consideration. 
The approach adopted by the MACHO collaboration was to 
divide their images into small sub-rasters and solve for the 
convolution kernel separately on each of these sub-rasters 
assuming no signific ant variations in the kernel occurred on 
their length scales (Alcock et al. 1999a). However, for the 
MOA images, we find small but significant spatial variations 
of the kernel even over small 500 x 500 sub-regions of the 
field. This prompted us to select the direct kernel modeling 
method which can easily model these spatial variations. 

Another consideration which influenced our choice of 
kernel solution method was that the PSF matching method 
requires high signal-to-noise empirical PSF models on both 
the reference and observation images. This would be difficult 
to achieve consistently using a small telescope such as ours 
at a site with only moderately good seeing. On the other 
hand, the direct kernel modeling method only requires the 
reference image to be of high signal-to-noise. 

An example of a subtracted image is shown in Fig. ^. 
The approach taken to subtract a MOA 2kx4k observation 
image from its corresponding reference image involved first 
selecting suitable stamps on the reference image. Typically 
~ 400 stamps were selected per 2kx4k field. An observation 
and reference image were then tiled into lkxlk sub-images 
which were small enough to allow the spatial variations in 
the kernel and differential sky background to be satisfacto- 
rily modeled using smooth functions. The solution for the 
convolution kernel and its associated spatial variations, to- 
gether with the differential sky backgrounds, was computed 
for each sub-region. The corresponding subtracted image 
was then formed using 



Ai(x, y) = i(x,y) - r * k(x,y) - b(x, y) 



(•») 



For a time series of observation images of a given field, a 
reference image was formed from the best seeing image or a 
combination of the best seeing images. The reference image 
is fixed over the entire time series. The stamps used for the 
kernel solution were also fixed for the time series of obser- 
vations. 



4 LA. Bond et al. 



A mask image mapping the distribution of saturated 
stars was also formed from the reference image. If a particu- 
lar star contained saturated pixels, the entire star profile was 
masked out. There is little useful data that can be obtained 
at the position of saturated stars. The mask was applied to 
all subtracted images to obtain final images. 

3.2 Data Analysis Environment 

In this work, two distinct modes of analysis were employed. 
The "online" analysis mode was used for real time reductions 
the large quantity of 2Kx4K images collected during obser- 
vations. A more detailed and rigorous "follow-up" analysis 
was then applied to individual events found in the online 
analysis. The online analysis system is shown schematically 
in Fig. |i[ The process involves the interplay of a large num- 
ber of image and numerical datasets. It is more helpful to 
regard the system as a data and image management sys- 
tem rather than a linear data reduction "pipeline" . The key 
features of the analysis system are as follows. 

3.2.1 Task Interface 

Each of the tasks linking the databases shown in Fig. ^ is 
controlled by a single Perl script. These scripts function as 
wrappers to lower level image and data analysis programs 
most of which are coded in C++. Each task can be run in- 
dividually or combined in other scripts to provide a higher 
layer of automation. This design provides considerable flexi- 
bility in operation, ranging from complete automation down 
to image-by-image interaction. 

3.2.2 Image Registration 

The image registration database provides information on 
the geometric alignment of observation images to their cor- 
responding astrometric reference images. The registration 
process involves cross-referencing resolved stars on the two 
images to achieve alignment. Each of the 2kx4k observa- 
tion images was sub-divided into lkxlk sub-frames. The 
brightest 1000 stars within each sub-region were found us- 
ing our implementation of the algorithm in the IRAF task 
DAOFIND. The star positions in each sub-frame were cross- 
referenced with star positions on the corresponding sub- 
frame of the astrometric reference using a smart algorithm 
which starts with the pre-supposition that the brighter stars 
on the observation image will be correspondingly bright on 
the astrometric reference image. A slower but very robust 
alternative algorithm was used as a backup if necessary. 
This combination of algorithms provided an extremely ro- 
bust frame-to-frame matching procedure that worked for all 
images taken under a variety of viewing conditions. For each 
observation image one finally obtained a list of 6000-8000 
cross referenced stars over the full 2kx4k field. The posi- 
tions of these stars on the observation and reference images 
were entered into the registration database. 

3.2.3 Image Alignment 

Geometrically aligned images were formed from flat-fielded 
images using information in the registration database. The 



alignment task was controlled by a Perl wrapper script for 
the IRAF tasks GEOMAP and GEOTRAN. This alignment 
process was used in the construction of subtraction reference 
images and in preparing images for the subtraction process 
itself. 



3.2.4 Reference Image Databases 

These databases contain the astrometric and image sub- 
traction reference images for each CCD of each observation 
field. The astrometric reference images define the coordinate 
system to which all observation images were geometrically 
aligned. The images used for both astrometry and subtrac- 
tion were selected from amongst the best seeing, highest 
signal to noise, and lowest airmass images obtained early in 
the 2000 Galactic Bulge observation season. We decided not 
to form subtraction reference images by combining several 
images. Instead we chose just one good seeing image per 
field as the reference. When combining large images such as 
ours, there is a danger of corrupting the linearity in the final 
image due to non-uniformities in the sky background across 
the chip. In image subtraction, it is crucial that the pixel 
values associated with the images are linear with respect to 
the incident flux.. 



3.2.5 From Raw to Subtracted Images 

The backbone of the image reduction flow involves pro- 
cesses that generate subtracted images from raw observation 
images. Intermediate image databases contain flat-fielded 
and aligned images. To relieve pressure on computer disk 
space, raw, flat-fielded, and aligned image databases were 
flushed periodically. However, all raw images were externally 
archived onto DLT tapes. When required, a given set of flat- 
fielded and aligned images could easily be regenerated using 
information in the registration database. On the other hand, 
all subtracted images were retained on disk as these were 
constantly being accessed by the photometry measurement 
processes. 



3.2.6 Identifying Objects on Subtracted Images 

The positions of variable stars on a given subtracted image 
exhibited positive or negative profiles depending on whether 
the flux had increased or decreased relative to the reference 
image. The subtracted images also contained spurious pro- 
files not associated with stellar variability. One type of spu- 
rious profile arose from the effects of differential refraction. 
Dependent upon the star's colour and the airmass at the 
time of observation, the position of its centroid was some- 
times displaced to a greater extent than other stars. The 
result was an imperfect subtraction at the star's position, 
giving a combined positive and negative profile. We did not 
attempt to correct for this effect, but instead filtered out 
these imperfect subtractions (see Alcock et al. 1999b for a 
discussion on how to deal with differential refraction). 

The first step in identifying objects associated with gen- 
uine variability was the application of the star finding algo- 
rithm described previously to the subtracted image to pro- 
duce an initial list of variable star positions. This list was 
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Figure 2. Road map of the on-line data analysis environment. The boxes represent the various databases. The dashed boxes datasets 
which are flushed periodically to relieve pressure on system disk space. The diamonds represent the analysis procedures connecting the 
databases. 



then screened by applying a second algorithm which exam- 
ined the profile at each position. Essentially, if the profile 
contained a large number of either positive definite or neg- 
ative definite pixel values, then the profile was deemed to 
be that of a genuine variable star. For each CCD on each 
field, newly detected objects were assigned running numbers 
according to discovery sequence. 

3.2.7 Photometry Databasing 

We did not maintain a database of light curves of all re- 
solved stars across the observation fields. Rather, we main- 
tained a database of variable star positions coupled with the 
database of subtracted images. Light curves were calculated 
only for objects of "current interest" such as newly iden- 
tified transient events. The variable star position database 
was updated periodically as new variables were identified 
and more refined positions were calculated. Since the sub- 
tracted images were retained on disk, a light curve at any 



given set of positions could be easily generated. In the online 
analysis, we extracted photometric measurements from the 
subtracted images using aperture photometry with an aper- 
ture radius of 6 pixels. In the follow-up analysis described 
in the next section, we used more a precise measurement 
technique. 

3.2.8 Event Selection 

We wished to detect new transient events online as opposed 
to periodic variables such as Miras, RR Lyrae stars, etc. A 
search for new events was generally carried out at the end 
of each night of good seeing conditions. For each field, the 
sequence of subtracted images from the previous night was 
combined to form a mean subtracted image which was exam- 
ined for positive and negative profiles due to stellar variabil- 
ity. The positions of these profiles were cross-referenced with 
those of variable stars previously detected to produce a list 
of "new" variables. The variable star database was then up- 
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dated by appending the new detections. Light curves for the 
new detections were generated by performing aperture pho- 
tometry at their respective positions on the entire sequence 
of subtracted images up to the current observation date. 
This practice of periodically carrying out new searches for 
events rendered the procedure sensitive to short timescale 
transient events. 

The light curves corresponding to all new detections 
were scrutinized by eye using a graphical tool. If a partic- 
ular light curve was clearly characteristic of periodic stellar 
variability, the associated position was flagged as hosting a 
"variable" star. On the other hand, if a light curve exhibited 
behaviour characteristic of a microlensing candidate or other 
transient event, it was flagged as "interesting". The interest- 
ing events were further scrutinized by examining both the 
subtracted and unsubtracted images at their respective po- 
sitions and checking for bad pixels, nearby saturated stars, 
bleeding columns etc. Interesting events were closely moni- 
tored throughout the viewing campaign. If any event could 
be clearly identified as a microlensing candidate or another 
tyoe of interesting transient event, a MO A Transient Alert 
was issued to the microlensing community and posted on 
the World Wide Web. | 

Early in the observation season, the number of light 
curves to inspect by eye was large but this number decreased 
rapidly as more and more periodic variables were flagged. 
As the observation season progressed, around 100-150 pre- 
viously undetected objects were found each time a search for 
new events was carried out. This number of corresponding 
light curves could easily be scanned by eye and classified. 
There is merit in being able to pick up events by eye in real 
time. No assumptions on the form of the variability of tran- 
sient events need be made in the selection process. As well 
as microlensing events, the possibility of detecting new and 
unusual transient phenomena arises with such a system. 



3.2.9 Computational Constraints 

The reduction of observational images using image subtrac- 
tion online is demanding on computer resources such as CPU 
time and disk space. The data analysis environment system 
described here was implemented on a dedicated 4-processor 
Sun Enterprise 450 computer located at the observatory site. 
Online disk space comprised 90 GB internal and 300 GB ex- 
ternal and was expanded from time-to-time. Raw images 
entered the analysis pool shown in Fig. |i| soon after acqui- 
sition and certification by the observing staff. A good night 
during the southern winter observing season yielded up to 6 
GB of raw images. In this case, the analysis took 24 hours to 
completely reduce the night's observations. However, it was 
not necessary to wait until all current observations were re- 
duced before searching for new events. 

Thus, despite the considerable demands on computer 
resources, we were successful in implementing an on-line im- 
age subtraction analysis system along with developing an 
effective system for issuing alerts of microlensing and other 
transient events in progress during 2000. 



t http:/ /www. phys.canterbury.ac.nz/~physib/alert/alert. html 




I 

Figure 3. Locations of the 14 MOA Galactic Bulge fields and 
positions of transient events detected in 2000 shown in Galactic 
coordinates. Candidate microlensing events detected during 2000 
are shown by the dots and the nova-like events are shown by the 
squares. The relative positions of the three CCDs are shown for 
field ngbl4. The MOA fields observed during 1998-99 are also 
shown. 



3.3 Selected Events and Follow-up Analysis 

A number of interesting transient events were detected in 
real time during the 2000 Galactic Bulge observation season. 
These are listed in Table [[]. Their positions are plotted in 
Galactic coordinates together with the MOA Bulge fields 
in Fig. H| The majority of the events listed are probably 
microlensing events, but the list also includes some nova- 
like events. Nine of the events listed in Table ^ were issued 
as transient alerts. 

The online analysis described in the previous sub- 
section involved a trade-off between analytical rigour and 
computing time. A detailed follow-up analysis was carried 
out on the events listed in Table nj. For a given event, 
400x400 pixel subframes were extracted from all flat- fielded 
observation images in the sequences. An astrometric refer- 
ence was again selected from amongst the best quality sub- 
frames and the geometric alignment was repeated. Since we 
were now dealing with small images, we formed the sub- 
traction reference by combining a selection of the very best 
seeing images. We required that the geometric alignment 
images selected in forming the reference involved small off- 
sets from the astrometric reference less than 5 pixels. This 
ensured any pixel defects are localised to small regions on 
the combined image after geometrically aligning the indi- 
vidual images. Furthermore, as with the online analysis, we 
imposed the strict requirement that all stamps on the refer- 
ence image used to solve for the kernel, be completely free 
of pixel defects. The image subtraction process was again 
carried out on all subframes. The use of small subframes 
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Table 1. Summary of All Detected Transient Events 





MOA TP) 




Coordinates (J2000.0) 




Alert Status 






Field 


PPFl 

\J ±J 


ID 


RA 


Dec 


MOA 




OGLE 




ngbl 


2 


ZUU I 


1 7- ^A- ^fi 


9Q-Q1 -A7 r. 
— zy.01.1 ( .0 


2000- 


-BLG- 


7 
1 








ngbl 


n 

z 


971 7 


1 7- ^7-fl7 Q1 




2000- 


-BLG- 


11 








ngbl 


q 
o 


797 
t z i 


1 7-^4-9Q 77 
J. I .oi.zy. 1 f 




2000- 


-BLG- 


3 








ngbl 


q 
o 




1 7- ^£-90 CM 
1 ( .oo.zu.y4: 


— zo.l 1 .^±0.0 














ngbl 


O 




J. ( .OO.UO.IO 
















ngb2 


2 




1 &-nn-i 9 ^fi 

lO.UU. 1Z.OU 


OQ.Q7-9Q Q 














ngb3 


2 


1 31 fi 


j-O.uo.uy.00 


— ou.ou.uu.o 


2000- 


-BLG- 


9 










i 


9snfi 

ZOUU 


1 7- c i r »-'} £ i 9(1 
1 1 .00.00. zu 


9C.1 n-1 7 1 


2000- 


-BLG- 










ngb4 


2 


91 Q7 
z ±y ( 


1 7- ^7-90 9fi 


— 97-4fv91 9 














n^t)4 


3 


159 


1 7-^7-zL7 90 
1 1 . J 1 .1 ( .zu 


— 97-^-^9 8 














ngb5 


i 


1 fi1 fi 


1 &-m -^n S7 

lO.Ui-.OU.O I 


90. c;Q.9fi Q 








2000- 


-BUL- 


ou 


ngb5 


i 


1 R9Q 


1 7- ^Q- ^8 1 1 
1 / .oy .00. 1 i. 
















ngb5 


]_ 


1 fifiS 
J-UUo 


1 k-fl1 -AA 7Q 

10.U1.t:4. 1 y 










2000- 


-BUL- 


io 


ngb5 


1^ 


1 R79 
i-U I z 


1 &-m -9fi 81 

lO.U-L.ZU.Oi- 


9R-^9-^4 7 








2000- 


-BUL- 




ngb5 


1 


1 fi73 
i-U f o 


18:01:06.71 


— 28:52:22.3 








2000- 


-BUL- 


■K) 




3 


1425 


18:03:54.78 


— 28:34:58.6 


2000- 


-BLG- 


12 


2000- 


-BUL- 


64 


ngb7 


3 


703 


18:10:55.62 


-29:03:54.2 


2000- 


-BLG- 


-8 








ngb9 


3 


841 


18:10:17.99 


-27:31:19.3 














ngblO 


1 


211 


18:08:07.97 


-26:13:14.4 














ngbll 


2 


1011 


18:11:51.18 


-26:26:48.6 














ngbll 


2 


1142 


18:11:28.31 


-26:15:05.8 


2000- 


-BLG- 


10 


2000- 


-BUL- 


-33 


ngbll 


2 


1146 


18:10:37.44 


-26:20:00.1 














ngbl2 


2 


1052 


18:14:47.42 


-25:32:53.6 














ngbl 3 


2 


1170 


18:16:56.78 


-23:29:50.9 


2000- 


-BLG- 


6 









made it easier to model spatial variations in the kernel and 
the differential background using smooth functions. 

On a given observation image, the profile on a sub- 
tracted image at the position of a variable star follows the 
point-spread function on the corresponding unsubtracted 
image. In the follow-up analysis the integrated flux over 
these profiles was measured using empirical PSF models. For 
each event, a high signal to noise PSF was constructed from 
the reference image by combining small sub-frames centred 
on a number of bright, isolated, and non-saturated stars near 
the event position. For each observation image, this reference 
PSF was convolved with the corresponding kernel model de- 
rived for this image and evaluated at the event position to 
obtain a normalized empirical PSF profile for the given event 
on the given observation image. This PSF was then rescaled 
to fit the observed flux profile on the subtracted image to 
provide a measurement of the flux difference. The associ- 
ated errors were determined in a self-consistent manner by 
applying the same measurement technique at the positions 
of constant stars in the field (which have flux differences of 
zero). The scatter in the flux difference measurements for 
these stars determined the frame-to-frame errors. 



3.4 Interpretation of Difference Imaging 
Photometry 

Difference imaging measures flux differences. This is a rela- 
tively unusual concept in astronomical photometry. To con- 
vert "delta flux" (Af 1 ) measurements onto a magnitude 
scale, it is necessary to determine a corresponding "total 
flux" . This is given by 



• CCD 1 




■ CCD 2 




• CCD 3 




-8 



-10 



-12 



-14 



MOA red 



F = F rcf + AF 



(6) 



Figure 4. Calibration of MOA custom red passband measure- 
ments to the standard I band magnitude scale for the three CCDs 
of the camera MOA-cam2. 



where F Ie f is the flux of the object on the subtraction refer- 
ence image. The reference flux may be measured by apply- 
ing a profile fitting program such as Dophot to the reference 
image. Alternatively, one can derive the reference flux by 
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fitting a theoretical profile, such as a microlensing profile, to 
the observed delta-flux light curve. 

There are some situations where it may not be possible 
to measure the reference flux, and hence the total flux, for 
a particular event. For example, the event may be strongly 
blended with a close bright star making a reference flux mea- 
surement impossible. In such a situation, classical profile 
fitting photometry will not be able to yield a reliable light 
curve but difference imaging photometry will at least yield 
a delta flux light curve with each measurement unaffected 
by blending. 

To obtain information in standard passbands, we used 
UBVI measurements of selected stars in Baade's Window 



provided by OGLE ( Paczyhski et al. 199£ ) that are con- 
tained in some of the MOA Bulge fields. The Dophot proce- 
dure was applied to the corresponding MOA reference im- 
ages and the resulting star positions were cross-referenced 
with those of stars with UBVI measurements. As seen in 
Fig. |i| we find that measurements taken in the non-standard 
MOA red passband correlate well with measurements taken 
in the standard I passband. If the total flux of a particular 
measurement obtained by image subtraction photometry is 
known, a corresponding standard I measurement can be de- 
termined by the following relations for the three CCDs of 
MOA-cam2 



f 25.8016 + 0.9704m rcd for CCD 1 
1 = 1 26.2964 + 0.9970m rcd for CCD 2 
1 26.9550 + 1.0316m rG d for CCD 3 



(7) 



where m rc d = — 2.51og(F 1 ) is the instrumental MOA red 
magnitude associated with the total flux measurement. 



4 EVENTS DETECTED DURING 2000 
4.1 Microlensing Events 

The observed flux difference for a microlensing event is given 
by 



AF = F base A(u(i))-F rcf 



(8) 



where Fbase is the baseline flux and F ie { is the flux on the 
subtraction reference image. In general the reference image 
contains some lensed flux and is not equal to the baseline 
flux. Even if the source is not visible at baseline, the best 
seeing images used to form the subtraction reference image 
may occur when the source is lensed. This is true for some 
of the microlensing events presented here. 

If the reference image was formed by combining iV fe f im- 
ages taken at times {ti}, the reference flux for a microlensing 
event is given by 



F„ 



-fbasc 

~N, 



— 5>(«(t0) 

rcf ^— ' 



(9) 



The amplification, A(u), expressed in terms of the dis- 
tance, u, of the lens star from the line-of-sight to the source 
s tar expressed in units of the Einstein radius is given by 
flPaczynski 1986|) 



A(u) = 



where 



+ 2 



Wii 2 + 4 



(10) 




*max 

Figure 7. Contour plot of \ 2 based on fitting a theoretical mi- 
crolensing profile in (A max , tg) parameter space for event ngb6- 
3-1425. The contours correspond to confidence levels of 68%, 95%, 
and 99.73%. 



u(t) 



t-t n 



t.E 



(11) 



Here u m i n is the minimum value of u during an event, £e 
is the characteristic Einstein radius crossing timescale, and 
tmax is the time of maximum amplification. The time £e is 
equal to te/vt where vt is the transverse velocity of the 
lens with respect to the line-of-sight to the source. 

We fitted the theoretical single lens light curve given 
by Eqns ^-|ll| to the data for candidate microlensing events 
using four parameters -Phase, Wmin, £e, and t ma x- The light 
curves are shown in Figs. H-n where the AF data points 
have been converted to amplifications using the fitted pa- 
rameters. We note that for some of the events, the baseline 
flux is below the observational threshold. However, we could 
still derive a best fit value for Fbase. This is essentially an 
extrapolation of the light curve from where the event is vis- 
ible down to the baseline level. In these events this value 
for .Fbase is used to calculate amplifications in plotting the 
light curves in Figs. ^-j^. We address the issue of measuring 
microlensing parameters for these types of events in detail 
below. 

The microlensing parameters £e and the peak amplifica- 
tion Amax (calculated directly from it m in using Eqn. Iicj) are 
the interesting physical quantities in a microlensing event. 
Both quantities can be measured for microlensing events 
with good sampling from the baseline to the peak. How- 
ever, for high magnification events sampled only around the 
times of peak amplification, it is possible only to measure the 
quantity t e g = te^min which is degenerate in the two param- 
eters of interest ( |Gould 1996 ) . We investigated the extent to 
which these parameters can be constrained by mapping x 2 
in the two dimensional parameter space corresponding to 
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Table 2. Fitted Parameters of Microlcnsing Candidates 

MOA ID ^4max £e / days t max I magnitude 

Field CCD ID lower best upper lower best upper JD— 2450000 baseline 
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Figure 8. Contour plot as in Fig. |?] for event ngb5-l— 1629. 



A m ax and £e. For each light curve, we fitted Eqns p|-|ll| for a 
sample of fixed (A max ,tB) pairs with the remaining param- 



eters, fbasc and t 

n 



allowed to float. 



In Fig. we show a contour map of % 2 in (/lmax,tE) 
space for the event ngb6-3-1425. The contours connect 
equal values of xLin + A* 2 for A x 2 = (2.3, 6.2, 11.8). These 
formally correspond to confidence limits of 68%, 95%, and 
99.73% respectively for the simult aneous estimation of two 
interesting parameters ( Avni 1976 ). Event ngb6-3-1425 can 
be described as a "classical" microlensing event where the 
source was resolved and visible at baseline and the light 



curve was well sampled both at its baseline and during its 
lensing phase. In this case both j4 ma x and £e are well con- 
strained. 

In contrast the parameter space of the event ngb5-l- 
1629 was not as well constrained as can be seen in Fig. ^. A 
close inspection of the baseline images for this event did not 
show any visible object at the position of the source while it 
was clearly visible during amplification. This event is clearly 
characteristic of one which is amplified from below the ob- 
servational threshold to above the threshold. Such events 
are only observable around the times of high magnification, 
and only lower limits can be inferred for /l m ax and tv- It is 
nevertheless noteworthy that one can infer that this event 
had a peak amplification well in excess of 70. 



Alcock et al. (200C) proposed classifying microlensing 
events found by DIA as either classical or " pixel lensing" 
events following the definition by Gould (1996) of pixel lens- 
ing events as microlensing of unresolved sources. A source 
may be unresolved if it is strongly blended with nearby stars 
or if it is too faint to be observed at baseline. Since DIA is 
unaffected by blending we suggest here a classification based 
just on amplification. Events would then fall into two types: 
those amplified from below the observation threshold and 
those amplified from above. For those events whose baselines 
are below the detection and hence measurement threshold, 
it is only possible to obtain lower limits on A mal and tv- We 
suggest that in DIA it is important to be able to distinguish 
between the two types of events as this affects the extent to 
which the microlensing parameters can be constrained. 

In Table |^ we summarize the results of this analysis for 
all 20 microlensing events showing best fit values for j4 max 
and £e along with 95% lower and upper limits where possi- 
ble. The fitted baselines fluxes were converted to I magni- 
tudes using the calibration information described in Section 
3.4. The parameters of one event, ngb3-2-1316, are poorly 
constrained as this event did not receive adequate temporal 
coverage. 
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Figure 9. Delta-flux light curve of MOA-2000-BLG-11 with 
standard and parallax microlcnsing fits. 



Table 3. MOA-2000-BLG-11: Standard and parallax microlens- 
ing fit ppraTTTBtCTs - 



Parameter Standard 



-^base 

to 

V 

e 

X 2 /dof 



14.16 

JD2451799.62 
43.04 days 
0.157 



18027.6/401 



Parallax 
14.78 

JD2451837.62 
69.69 days 
0.224 

42.52 km s" 1 
-76? 20 
1071.6/399 



Of the 20 events listed in Table [2j, half of them have 
best fit peak amplifications around 10 or more, and six have 
lower limits of A max greater than 10. This is a large fraction 
of high magnification events when compared with the cor- 
responding fractions from the MACHO and OGLE surveys. 
Since these groups use profile fitting photometry in their 
on-line analysis, their measurements are in general affected 
by blending and therefore the derived values of 4 maj . tend 
to be underestimated. More importantly, DIA is sensitive to 
microlensing events whose source stars are too faint to be 
detected at baseline. 

The large fraction of detected events with high mag- 
nification has useful implications for studies of extra-solar 
planets, low mass stars and stellar atmospheres, as noted in 
the introduction. 



4.2 MOA-2000-BLG-11: A Microlensing Parallax 
Event 

The light curve of event ngbl-2-2717, which was alerted as 
MOA-2000-BLG-11, is not well described by the standard 
microlensing profile given by Eqns. 7-8. The light curve ex- 
hibits an asymmetric profile which is characteristic of paral- 
lax microlensing. Here, a distortion arises from the circularly 
accelerating motion of the Earth around the Sun. When in- 
cluding this motion, the distance between the lens star and 



Table 4. Summary of the characteristics of nova-like events de- 
tected in 2000. The column, A7, denotes the change in I magni- 
tude during the rise of the event from its quiescent state to its 
high state. 





MOA ID 




peak I 


Al 


h/2 


Field 


CD 


ID 






days 


ngblO 


1 


211 


16.64 


>2.4 


5 


ngbll 


2 


1011 


15.78 


>3.3 


3 


ngbll 


2 


1146 


16.75 


>2.3 


7 


ngbl3 


2 


1211 


13.95 


>5 


4 



the source star (given by Eqn. 8 for the constant velocity 
case) contains two extra parameters: the angle, 0, between 
the transverse motion of the lens and the transverse speed, 
v, projected onto the solar position. Th e full mathe matica l 
details can be foun d in Alcock et al. fll995| ), Mao ( 1999 ), 
and Dominik ( |l998; ). 

The light curve of MOA-2000-BLG-11 is shown in 
Fig. |^ along with the best fit standard constant velocity 
and parallax microlensing models. The derived parameters 
are given in Table [| At 75.1 days, the duration of MOA- 
2000-BLG-ll is rather short compared with the durations 
for previously detected parallax events of 111.3 d ays for the 
first parallax event found by the MACHO group (Alcock et 



al. 1995[ ), 118.1 days for OGLE-99- CAR-1 (|Mao 1999[), a nd 



156.4 days for OGLE-00-BUL-43 flSoszyfiski et al. 200l| ). 

Measurements of microlensing parallax events partially 
resolve the degeneracy in the mass, lens distance, and trans- 
verse velocity contained in the standard microlensing equa- 
tion. Assuming negligible velocity dispersion in the disk and 
the bulge, the transverse speed projected onto the solar po- 
sition can be written as 

V = ^Lkms- (12) 
1 — x 

where x — Di/D s with Di being the distance to the lens 
and D a t he distance to the source star. Following Alcock et 
al. (1995) the mass of the lens, Mi, can written as 



Mi 



l~x v 2 tjc 2 
x 4GD S ' 



(13) 



Substituting the parameters obtained from our parallax fit 
and assuming D s = 8 kpc we obtain Di ~ 1.4 kpc and 
Mi ~ 0.3Mq for the lens star. It would then contribute very 
little light through the course of the event. It should be noted 
that, provided the lens star is constant, the contribution 
of any flux from the lens star is canceled out in delta-flux 
measurements. 



4.3 Nova-like Events 

Some events with light curves showing nova-like profiles were 
found. We present them in Fig. [jlj and Table ^ to illustrate 
the potential for detecting new types of objects. None of the 
events was visible in the MOA images during their times 
of minimum light. Event ngbl3-2-1170 (alerted as MOA- 
2000-BLG-6) is clearly characteristic of a dwarf nova. This 
event rose by more than five magnitudes in less than three 
days peaking at I — 13.95. Event ngbll-2-1146 also ap- 
pears to be the tip of a dwarf nova decay just above the 
detection limit. Event ngbl0-l-211 is also nova-like but the 
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Figure 10. Light curves of nova-like events detected during 2000. 



decay is not as smooth as the other events. At the posi- 
tion of event ngbll-2-1011, the MACHO database contains 
an object whose corresponding light curve undergoes flaring 
episodes about once per year (D. Bennett, private commu- 
nication) . 



4.4 Asteroids 

As well as measurements on variable stars, difference imag- 
ing analysis can be employed to study moving objects. In 
Fig. [H] we show subimages of three exposures of MOA field 
ngbl on CCD number 1 taking during one night. The subim- 
ages include the original unsubtracted images along with 
their corresponding subtracted images. The exposures were 
taken over a time interval of approximately 1.6 hours. An ob- 
ject moving across the field of view can clearly be seen. The 
positions and magnitudes of this object measured on each 
image are listed in Table ||. A crosscheck with Guide 7.0 
revealed asteroid (3680) Sasha at these positions and times. 

The associated transverse velocity across the sky for 
this object is 19.2" hour -1 . Frequent sampling of the fields 
are required if one is to track the motion of such fast moving 
objects during a night. Sampling the fields once per night 
would merely result in spike events at asteroid positions. 
The MOA data base now contains subtracted images derived 
from a large number of nights where the sampling rate was 
three or more times per night. We are currently exploring 
ways of exploiting this data base to study asteroids. 



Table 5. Astrometry and Photometry of Asteroid 



JD-2450000 

1675.066760 
1675.133438 
1675.197859 



Coordinates (J2000.0) 
RA DEC 



17:56:23.33 
17:56:21.28 
17:56:19.28 



-29:54:51.4 
-29:55:06.5 
-29:55:21.0 



I magnitude 

15.33 
15.19 
15.28 



5 SUMMARY 

It has been demonstrated that real-time difference imag- 
ing of an extensive database is possible. The procedure em- 
ployed by MOA is a database management system rather 
than a linear reduction pipeline. A variety of levels of au- 
tomation are possible but the final selection of events is by 
human intervention. 

To minimize hardware constraints, two reduction pro- 
cedures were developed, one for online analysis of large num- 
bers of large images, and one for follow-up, accurate analyses 
of selected events. As our on-site computing power expands, 
we plan to transfer elements of our follow-up procedure to 
the online analysis. We also plan to increase the conversion 
rate of internally monitored candidate events to actual alerts 
issued to the microlensing community. 

We have demonstrated that a high fraction of events de- 
tected by difference imaging have high magnification. This 
should have useful implications for studies of extra-solar 
planets and faint stars provided the peaks of future events 
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Figure 11. Superposition of sub-regions of 3 subtracted images 
taken on one night showing an asteroid trail. 



are monitored photometrically and/or spectroscopically by 
the microlensing and other astronomical communities. 

Our observational strategy has been to monitor a mod- 
erate area of the Galactic Bulge, 17deg 2 , fairly frequently, 
i.e. up to six times per night. This is a compromise between 
survey area and sampling rate that we plan to continue. It 
has enabled the detection of some unusual types of events, 
including nova-like transients and an asteroid, as well as mi- 
crolensing events of high magnification. 
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